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FOREWORD 


This  Geophysics  Research  Paper  is  published  in 
three  separate  volumes. 

In  Volume  I,  the  theory  of  large-scale  atmospheric 
diffusion  is  developed.  The  application  of  this  theory  to  air 
trajectories  is  presented  in  Volumes  II  and  III, 

Volume  II  contains  probability  tables  for  various 
constant  values  of  mean  zonal  wind. 

Volume  III  consists  of  probability  tables  for  North 
America  and  Eurasia, 
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ABSTRA  CT 


G.  I.  Taylor's  theory  of  diffusion  by  continuous  movement  is 
adapted  to  motions  on  the  scale  of  the  general  circulation.  The  re¬ 
sultant  theory  pertains  to  diffusion,  by  large-scale  eddies,  of  air 
particles  constrained  to  move  on  a  constant-level  surface.  This 
theory  provides  a  means  for  determining  the  probability  field  as  a 
function  of  time  for  a  particle  originating  from  a  given  point  on  the 
surface  of  a  sphere. 

The  bivariate  normal  density  function  describes  the  prob¬ 
ability  of  a  displacement  with  components  x  (West  -  East)  and  y 
(North  -  South).  This  function  reduces  to  the  circular  normal  form 
when  suitable  empirical  and  theoretical  simplifications  are  intro¬ 
duced,  concerning  the  mean  ?.onal  wind,  the  standard  deviation  of 
displacement  components,  and  the  correlation  between  displacement 
compoiuMU  deviations. 

This  dei»aity  function  is  integrated  in  polar  form  and  mapped 
on  a  projection  of  tlie  earth's  surface  by  means  of  a  suitable  spher¬ 
ical  transfur matiotv.  The  resulting  form  of  the  distribution  function 
is  applicable  to  problems  of  atomic  fallout. 

In  applying  diffusion  theory  to  balloon  operations,  llte  concept 
of  downstream  probability  density  function  is  introduced.  Thi.«  func¬ 
tion  defines  the  probability  density  with  respect  to  latitude  of  an  East- 
West  displacement  at  least  as  large  as  a  specifivul  value,  occurring 
within  T  (Lays,  ‘rabies  of  the  downstream  probability  function  are 
presented  in  Volumes  II  and  lU  for  (a)  various  cvmslani  values  of 
mean  wind  and  (b)  Xortli  .America  and  Eurasia  for  January,  April, 
July,  and  OcUiber  at  <id,00l)  to  80,000  feet. 
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THEORY  OF  LARGE-SCALE  ATMOSPHERIC  DIFFUSION 
AND  ITS  APPLICATION  TO  AIR  TRAJECTORIES 


1,  INTRODUCTION 

In  recent  years  the  constant-level  balloon  has  developed  into  a  vehicle 
of  xTvajor  scientific  importance.  Its  payload  may  consist  of  a  variety  of  jxieas- 
uring  devices  which  utiliae  appropriate  communications  systems  to  telemeter 
information  back  to  earth.  In  many  cases  this  kind  of  data  can  be  obtained  in 
no  other  way.  Indeed,  even  the  patl'S  of  these  balloons  are  of  great  scientific 
interest,  since  they  represent  solutions  to  the  equations  of  large-scale  atmos¬ 
pheric  motions  m  a  Lagrangean  system  of  coordinates. 

Une  of  the  most  attractive  features  of  the  constant-level  balloon  is  a 
direct  coiusequence  of  its  aiKe.  U  is  large  enough  to  be  inseiusitivy  to  small- 
scale  eddies,  while  responding  freely  to  the  large-scale  motions  of  the  atmos¬ 
phere,  Parado.HicaUy,  it  is  tins  very  freedom  of  movement  which  imposes  the 
most  severe  operational  limitatitui  on  the  hallomi  .system,  namely  uncontrol - 
lability.  The  horizontal  path  can  be  coturoUed  to  a  small  extent  in  two  vvays. 
(a)  by  varying  the  time  of  launching,  anti  thus  predetermining  the  initial  hori¬ 
zontal  Velocity  at  altitude,  d  fh)  by  v.,i»  ying  the  altitude  during  flight.  Except 
tor  these  measures,  knowledge  of  the  future  positions  of  balloons  must  largely 
depend  upon  meteorological  prediction. 

Experience  has  sh«twn  that  for  lime  periods  great er  (h.in  a  few  tiays  con¬ 
ventional  pri'dicliim  methods  are  unreliable,  r'or  these  longer  time  intervals, 
statistical  metlo.tds  must  be  used.  I'o  be  valid,  sueli  rnct'nods  must  be  based 
on  a  set  of  fundamental  sl-itisties  of  large-scale  atmospheric  motions.  Tiiis 
IS  hut  one  example  of  the  need  for  .sucli  statistics. 

NUTlf;  Author's  m.unsscripl  approved  Id  March  l’?58. 


li.i  vv'  ijri)po;;i'cl  lor  oonibitu’d  lu-l.vorks  ot'  b;illoons, 

l.i’Moh  ;  I'roni  ,so\  eral  ^ilo;:  .it  i)ri.'iK'U'rini.'U'cl  tiiao  intiTval.s.  Suoh  oN,pi'ri- 

I. '.-  as  .iro  noi'i'saarily  ooiUly  aad  must  bo  thoruu;’,hly  ovaluaU-d.  To  porlorm 
..a  lilt  ot  ■  ip  out  job  ot  ovaluatioa,  ouo  must  apaiu  liavt'  aoci'ss  to  a  relbablo  set 

oi  at. Ill  iLos  on  atmusnhoric  tm)lions.  lia  I'aot,  th.-  lack  ot  such  intormation  has 
soriou  y  dotorrod  reali^atiui^  of  tho  full  potential  value  of  b.illoon  systenas. 

Tiis  paper  presents  a  sot  of  statistics  on  atmv)splieric  trajectories  ifi 
a  lorm  vhich  has  been  found  most  useful  for  solvinp  the  problems  discussed 
above  .1  id  other  related  problems. 

the  most  straightforward  approach  to  the  problem  of  distribution  of 
trajectory  end-points  is  to  compute  ti  large  number  of  trajectories,  having  a 
mumn  a  origin,  from  a  series  of  upper-level  chtirts.  Here,  either  the  ob- 
s  oo  vinds  or  the  computed  geostrophic  winds,  or  both,  are  used  as  basic 
liat  i.  1  he  accuracy  of  such  computations  will  depend  to  some  extent  on  the 
lime  ii.  erval  between  successive  mtips.  A  numb>  r  of  such  studies  have  been 
n  ido.  nd  these  are  useful  in  giving  at  least  a  gros;>  estimate  of  the  diutribu- 
tion  fur.,  lions.  However,  itt  addition  to  the  magnitude  <.*t  this  task,  Ihgre  are 
tier  rious  limitations  to  this  a^jproach; 

a.  U  is  exiretiu'ly  difltcull  to  ctmtpuie  aceurale  tra- 
,■  cturies. 

b.  With  the  type  ol  data  at  our  dtssjmsal,  it  is  almost 
;mp<jsstble  to  Achieve  a  .sample  sufficiettily  large  to  be  cun- 
■  dere<i  representative. 

Keen  m  iho.se  few  c-ises  where  the  m'Mt  iwa  Hmtta' 
t.  -ns  CiUt  be  overeome,  the  appHcabiHly  of  the  results  is  re- 
''rutled  to  A  particular  region  in  space,  a  particular  height, 

•-  id  a  pa rticuiar  season, 

.V:th  regard  to  the  first  limUation.  it  is  only  fair  to  stale  that  a  sufficiently 
large  o .1, mpl e  would  overcame  this  restriction  if  there  wore  no  systematic  com- 
p-Uaii.-.  errors. 

iiie  .Moby  Hick  flights  are  an  important  source  of  data.  Although  many 

J,  Ur.. ire  s  oi  balloons  have  been  flown,  only  a  smalt  percentage  of  these  flights 
•ts  e  r  uffu'ienlly  long  duration  to  be  useful.  Thus  the  la.sl  two  limitations 

r<  nain  in  force,  Kurtheinirire,  there  is  a  definite  hi.vs  in  these  data  bc- 
t ue  launchings  were  maivdy  confined  to  selected  situations  when  the  wcat- 
rHvt  .ere  strongest. 


From  the  foregoing  discussion  it  is  evident  tliat  an  adeciuate  theory 
of  diffusion  on  the  scale  of  the  general  circulation  is  needl'd.  This  paper 
attempts  to  derive  such  a  theory.  This  is  by  no  means  a  complete  theory, 
since  certain  simplifying  assumptions  and  approximations  have  been  in¬ 
troduced,  The  criterion  adopted  in  introducing  these  approximations  was 
to  maintain  compatibility  between  the  degree  of  refinement  of  the  theory 
and  the  quality  of  the  input  data. 

Any  theory  of  diffusion  of  air  particles  by  the  general  circulation  ed¬ 
dies  must  provide  a  means  for  determining  the  probability  field  as  a  function 
of  time  for  a  particle  released  from  a  "iven  point  on  a  sphere's  surface. 
Initially,  the  probability  of  the  particle  s  location  is  100  percent  at  the  ori¬ 
gin  and  0  elsewhere.  As  time  passes,  the  particle's  location  is  specified 
by  a  probability  field.  After  a  very  long  period  ot  time,  the  probability 
field  tends  toward  umforinity  over  the  sphere's  entire  surface  and  the  par¬ 
ticle's  origin  ceases  to  be  important. 

G.  I,  Taylor's  theory  of  diffusion  by  continuous  movement^  is  a  solu¬ 
tion  of  this  problem  for  small-scale  eddies.  The  theory  presented  here  is 
an  adaptation  of  I'aylor's  theory  for  motions  on  the  scale  of  tlie  general  cir¬ 
culation,  in  developing  our  theory,  we  incorporated  results  obtained  from 
empirical  studies  of  the  trajectory  distribution  function. 

The  probability  density  function  can,  of  course,  be  integrated  witli  re¬ 
spect  to  .space  and  time  in  a  number  of  ways  to  obtain  the  desiri'd  distribution 
function,  in  planning  tlie  p c  c,>j entatiou  of  tlieorelic.il  cal I'ulations  ,  we  found 
that  each  specific  class  of  application  required  a  particular  metliod  of  inte¬ 
gration  to  obtain  the  most  useUil  form  of  the  distribution  function.  One  pre¬ 
sentation  is  included  'vvliica,  while  not  very  useful  for  con.sUilU-level  balloou 
operations,  may  have  imporlanl  applications  ti.i  ilie  problem  of  tUffusion  v)f 
paniculate  matter,  sucli  as  fallout  from  atomic  weapons. 

d.  l-MAOHAllli,lTY  'i'lil'.'UKV  Ol-'  GhS  I  *1  .A  G  f,  M  l-.'iS  I'.S 

Gonsider  an  air  iMiiicte  conslniined  lo  move  on  .i  c  on  stvi  iil  -  I  ev  el  surlace. 
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This  particle  is  located  at  0  at  time 
t  =  0  (Figure  1).  What  is  the  proba¬ 
bility  of  the  particle  being  located  at 
point  P,  at  t  =  T?  (Note  that  in  this 
statement  of  the  problem  we  prescribe 
only  the  end-point  of  a  displacement, 
not  the  path.)  Making  the  usual  assump¬ 
tion  that  the  two  displacement  compo¬ 
nents,  X  and  y,  are  normally  distrib¬ 
uted  about  their  respective  means,  X 
and  Y,  the  probability  of  displacement, 
OP,  may  be  expressed  by  the  equation 
for  the  bivariate  normal  density  func¬ 
tion:  i" 


P 


!  y 


X 


Figure  1.  Displacement  of  an  air 
particle 
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0  d 
X  y 


\  d 


(i) 


where  and  6^.  are  the  standard  deviation  of  x  and  y,  respectively,  and  p  is 


the  coefficient  of  correlation  between  the  quantities  x’ 
liy  definition 


X  = 


T 

f  U  dt 


:2) 


and 


Y 


(x”X)  and  y'  -  (y-Y). 

T 

I  V  dt 


(i) 


where  1)  and  V  are  components  of  the  dimatological  mean  wind  encountered  by 
the  particle  in  its  pnili. 

In  the  stratosphere  and  upper  troposphere,  V  is  small  and  difficult  lu  dv'- 
tormine.  This  is  especially  true  in  tnidille  lutllvuies.  In  fact  it  may  be  assumed, 
with  little  loss  of  gener.ilily,  that  the  mean  wind  in  thc.se  regions  is  /.ojtal,  l.et 
us  flow  more  closely  examine  the  mean  wiiul  variatiuu. 

Tahle  1  is  a  typical  satuple  of  mean  upper-U'vel  v.  inils,  taken  i  rum  thv 
Haiuibuuk  vif  (leoohysics.  “  Note  lliat  the  vanalion  is  coiusulerably  gre.iier  wuii 
lo-spect  to  latitude  than  with  respect  to  longitude. 
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Table  1;  Mean  wind  speed  (knots),  50,000  ft.,  January 


LATITUDE 

WEST  LONGITUDE 

LONGITUDE  RANGE 

120 

100 

80 

60 

70 

34 

33 

33 

34 

i 

60 

36 

29 

26 

29 

10 

j 

50 

38 

41 

44 

43 

6 

40 

71 

84 

83 

67 

17 

30 

61 

58 

48 

37 

24 

LAT  KANGE 

37 

55 

57 

38 

It  is  evident  that  if  we  confine  our  discussion  to  specific  geographical  regions, 
to  a  reasonable  approximation,  we  need  only  consider  the  latitudinal  variation 
of  the  mean  wind.  Therefore,  we  introduce  the  following  modeling  approxima¬ 
tions; 

Y  =  V  =  0  (4)  and  U  =  U._  (())  .(5) 

where  i,  j,  and  k  are  height,  season,  and  geographical  region  (such  as  North 
America,  Eurasia,  etc.  ),  respectively,  and  0  is  latitude. 

It  can  be  shown  theoretically  that  p(x',  y'),  the  coefficient  of  correlation 
between  the  displacement  components,  is  proportional  to  p(u',  v'),  the  corre- 
uition  coefficient  between  the  wind  components,  where  u*  =  (u-U).  and  v'  s  v. 
Brooks,  ft  al,  ^Buch,^  and  Couid.  conducted  independent  studies  of  statist¬ 
ical  pre;  eiTies  of  the  winds  in  the  free?  atmosphere  over  broad  regions.  All 
three  authorities  agv?u>  that  p{u',  v')  is  small.  Table  i  is  u  typical  example 
extracted  from  Court's  paper. 

Table  Z:  tloefticiont  of  correlation  between  u  and  v 


HT. 

(km) 

PA  N A  MA 
Winter 

CANAL  ZONE 
Summer 

MT.  CL 
Winter 

EMENS.  MICH. 
Summer 

GOOSE  BAY,  l.AB. 
Winter  Summer 

ft 

.  000 

.  124 

.054 

.  05  1 

-.119  .022 

It) 

.  075 

.  04  1 

.091 

.  050 

-.087  .02? 

20 

.059 

174 

.  129 

154 

.11.1  -.  063  1 

In  Eq.  (1),  P  ,  which  appears  in  both  denominator  and  exponent,  must 

be  compared  to  unity.  Since  the  first  and  last  terms  in  the  exponent  within 

the  brackets  are  of  the  same  order,  the  value  of  p  within  the  bracket  must 

be  compared  to  unity.  From  a  complete  study  of  the  evidence  it  becomes  ap- 
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parent  that,  compared  to  1,  p  is  small  and  p  is  very  small.  Thus  we 
introduce  the  assumption  that 

P  -  0  (6) 


The  time  dependence  of  and  0^  is  fully  discussed,  in  Appendix  III 
where  the  equations  ^y  ~  developed.  It  is  also 

shown  in  Appendix  III  that  empirical  evidence  warrants  the  assumption  that 


G  (t)  =  G  (t).  Therefore,  it  is  evident  that  6  =  ciO  where  q  is  the  ratio 

X  V  y  ‘  X  ^ 

6v 


"u 


Consider  now  the  quantities  6^^  and  G  .  With  the  sample  sizes  gener¬ 


ally  available,  G^  can  be  determined  with  much  greater  precision  than 


'I’hereiore,  it  is  reasonable  to  assume  that  any  error  in  the  estimated  value 


of  q  will  be  reflected  in  G^  rather  thai^  in  d^.  Let  us  then  investigate  the 


consequences  of  assuming  that  q  =  i.  If  q  is  actually  greater  than  1,  then 
the  dispersion  in  tlie  y-directioit  is  greater  than  we  liave  assumed.  On  the 
other  hand,  if  q  Is  actually  smaller  than  1,  then  we  have  overestimated 
tlte  dispersion  in  the  y-direction. 

In  applying  tlie  theory  to  balloon  systems,  we  are  usually  interested 
in  obtaining  the  greatest  possible  concentration  (tliat  is.  least  dispersion) 
ot  displacements.  .As  will  be  shown  later,  the  results  of  several  studie^s  in- 
dii-ate  tlial,  in  geiu'ral,  q-  1.  l  lms  it  may  be  seen  that  any  errors  resulting 
from  the  q  "  1  assumption  will  be  in  tlie  eonservalive  direcluni  when  the 
llieory  is  applied  to  balloon  systems.  W^^  will  now  discuss  llie  magmluile  ot 
iheSt'  errors. 

i.d 


extensivelv  inv estiiia t lul  G  aiul  u  and  cuiuduiied 

•  ••  VI  V 


L  rook  s ,  I'l  a  I 

that  there  was  no  signifiinnl  dlffereiu'e  bolwetui  ll  e.se  ouautities  in  tlie  free 
a imos (die r e .  lUich  conducted  .i  detaiU'd  iuvesligaiiun ,  on  .a  henuspherii.: 
b.i,si,s.  oi  the  wind  statislu  s  ft.ir  a  single  year,  llis  findings  agree  witli  thosv' 


of  Brooks'  in  middle  and  high  latitudes.  However,  he  found  that  in  low  lati¬ 
tudes  q  is  significantly  less  than  1,  .  Court's^  more  recent  studies  show  that 
there  is  a  consistent  variation  of  q  with  elevation  at  widely  separated  sta¬ 
tions.  The  ratio,  q  is  1  at  the  Jet  Stream  level  (about  40,  000  ft)  and  de¬ 
creases  gradually  with  increasing  elevation.  Court  also  found  that  this  effect 
is  very  pronounced  at  low  latitudes,  but  it  is  much  smaller  at  middle  and 
high  latitudes.  Let  us  confine  our  attention  to  middle  and  high  latitudes  and 
to  elevations  between  40,  000  and  80,  000  feet.  All  authorities  menticr.cd  above 
agreed  that  the  error  in  assuming  6^  =  6^  is  negligible  at.40,  000  feet.  Our 
own  studies  show  no  significant  difference  between  6  and  6  ,  even  at  80.  000 
feet.  These  results  are  not  necessarily  in  conflict  with  Court's  data.  He 
used  data  from  single  stations,  whereas  our  estimates  were  based  on  popu¬ 
lations  of  winds  taken  along  singl platitudes  but  from  a  number  of  adjacent 
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longitudes.  By  definition,  =  v  -  (v)  ,  It  has  already  been  shown  that 

the  regional  va'ue  of  (v)  is  0  (Eq.  4  ).  For  a  single  point,  however,  (v)^ 
is  not  negligible,'  Therefore,  other  things  being  equal,  the  value  of  d^  will 
always  be  larger  for  a  region  than  for  a  single  point  within  the  region. 

The  foregoing  discussion  leads  to  the  assumption  that 

6y=d^  =  d  .  (7) 

VVe  may  now  rewrite  Eq.  (i)  in  various  alternate  forms: 


P  (x,  y,  t)  =  p  (X,  t)  p(y,  t)  = 


- - e.xp 

2’-  d“ 


/  t  \  ^  i  ^ 

(x*)  4  y 
?.d^ 


p(x,t)  =  - 

V27T  d 


P(y.t)-  T— -  exp  --^3  .  (10) 

\/27Td  [  2d^ 

Wi>  may  write  Eq.  (vS)  in  polar  form  by  introducing  -  (x')*^  +  y*^  and 
ialtigraling  from  0  to  2'T  ; 


The  preceding  equations  are  familiar  forms  of  the  circular  normal  density 
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function.  Brooks,  et  al,  ’  used  these  equations  to  compute  the  probability 
of  winds. 


3.  ANALOGY  BETWEEN  THE  WIND  PROBLEM  AND  THE  DISPLACEMENT 
PROBLEM 


The  displacement  problem  can  be  viewed  as  a  logical  extension  of  the 
wind  problem.  We  will  compare  the  solutions  to  these  two  problems  at  every 
step  to  illustrate  their  relationship.  In  this  analogy  we  will  retain  the  sim¬ 
plifying  assumptions  and  approximations  which  have  already  been  established. 

The  appropriate  wind  equations  which  correspond  to  Eqs.  (8),  '  (10), 

and  (11)  are: 


p(u,  v)  =  p(u)  p(v)  = 


ZlTd.' 


exp 


-  -A 

(u')^  +  v^ 


2d.- 


■  • 


(8a) 


P(u)  = 


1 


.  exp 

V2ir  0. 


M. 

26.^ 


(9a) 


p(v)  =  z-h  exp 

V  ^  n'ci. 

aiui 

p(c)  =  y  ,  exp  j  -  j 

o,‘-  2d."  I 

I  III 


(iOa) 


(lla) 


where  d.  is  the  standard  wind  deviation  and  =  (u')*^'  +  . 

We  note  that  the  indopendur.t  variables  fur  the  wind  equations  are  u  and 
v;  but  for  the  displacement  equations,  in  addition  to  x  and  y,  the  variable,  t, 
also  appears.  In  general  botli  the  wind  and  displacement  probability  density 
funetioi\s  are  completely  specified  when  two  primary  paranieters,  U  and  6.  , 
are  given.  However,  the  empirically  determined  relationsliip  between  U  and 
d.,  devi'loptui  in  .‘\ppendix  11,  reduees  llu'  required  number  of  primary  para- 
tnelers  to  a  single  one,  natmdy,  U.  We  define  secondary  parameters  as  those 
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variables  which  represent  the  time  and  space  dependence  of  the  mean  wind  it¬ 
self.  In  this  paper  the  secondary  parameters  are  height,  season,  region  (in 
the  sense  defined  previously),  and  latitude  for  the  wind  statistic  and  similarly, 
height,  season,  region,  and  initial  latitude  for  the  displacement  statistic. 

In  the  wind  problem  the  specified  wind  (u,  v)  can  be  represented  as  a  two- 
dimensional  vector  in  plane  surface.  However,  the  vector  x,  y  must  be  mapped 
on  the  earth's  surface  which  is  spherical.  Therefore,  at  some  point  in  the  so¬ 
lution  of  the  displacement  problem  we  must  introduce  a  transformation  to  spher¬ 
ical  coordinates.  One  additional  complication  arises  which  is  confined  exclus¬ 
ively  to  the  displacement  problem.  Since  the  component  y  represents  a  change 
in  latitude,  and  since  both  U  and  are  functions  of  latitude,  we  must  somehow 
incorporate  in  the  probability  function  a  continuously  varying  population  of  winds. 
Just  how  we  overcome  these  difficulties  will  be  discussed  in  the  next  section. 

4.  POLAR  INTEGRATION  OF  THE  PROBABILITY  DENSITY  FUNCTION 

Any  desired  form  of  the  distribution  function  may  be  obtained  by  an  ap¬ 
propriate  integration  of  Eqs.  (8)  to  (11)  and  (8a)  to  (11a).  A  useful  form  of  the 

wind  distribution  function  is  derived  from  integration  of  Eq.  (Ha).  A  rather 
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complete  discussion,  of  this  function  will  be  found  in  Brooks,  et  al.  * 

Integrating  Eq.  (Ha)  with  respect  to  c: 


This  function  may  be  interpreted  as  follows:  Any  wind  with  probability,  P,  is 
the  vector  sum  of  U,  the  vector  mean  wind,  and  a  deviation  vector  of  magni¬ 
tude,  C,  where 

tn  ^  (13a) 

^1  -  P(C)  ^ 

An  analogous  polar  integration  of  the  displacement  density  fuuction(Eq,  11) 
is  useful  in  certain  types  of  problems.  To  demonstrate  the  method  of  computing 
this  form  of  distribution  function,  we  will  apply  the  theory  to  a  specific  example. 
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sr^ 

,v'>' 

,v\%’ 
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m 


The  problem  we  have  selected  is  the  mapping  of  the  0.  50  probability  curve 
for  the  following  specifications: 

Region;  North  America 
Height;  50,  000  Feet 
Season:  January 

Initial  Latitude:  TO'^N 
Time:  1  Day 

The  profiles  of  U  and  6^  are  shown  in  Figure  2.  As  a  preliminary  to  the 
discussion  of  the  displacement  function,  we  solve  the  0.50  probability  prob¬ 
lem  for  the  wind  occurring  at  40^N.  We  note  that  the  mean  wind  at  40°N  is 
75  knots  and  the  standard  deviation  is  35.  8  knots  (Figure  2).  From  Eq.  (13a); 

C  =  35.  8  [  In  4  J  =  35.  8  x  1.  176  =  42  kts.  (14a) 

The  solution  is  illustrated  in  Figure  3.  The  probability  of  occurrence  is  at 
least  0.50  for  all  winds  which  can  be  represented  by  vectors  with  a  common 
origin  at  0,  0  and  whose  end-points  lie  within  the  designated  circle.  Thus 
the  circle  of  radius,  42  knots,  and  origin  at  the  end-point  of  the  vector,  U=75 
knots,  is  the  locus  of  the  limiting  vector  end-points  for  a  probability -of  0.50, 
The  solution  consists  of  two  independent  parts;  (a)  a  mean  wind  vector,  U, 
and  (b)  a  deviation  wind  vector,  u',  cf  magnitude,  C.  Figure  3  also  clearly 
shows  that  the  component  u  =  (U  +  u'). 

For  the  displacement  distribution  function  we  seek  a  sitnilar  solutioji, 
but  we  further  require  that  the  0.50  probability  curve  be  mapped  on  a  spher¬ 
ical  earth.  For  this  purpose  i-  is  convenient  to  use  one  degree  of  latitude  on 
liiu  earlii's  surface  as  a  unit  of  length,  and  one  day  as  a  time  unit.  U  is  also 
convenient  to  adopt  a  coordinate  system  in  which  /  is  the  longitude  coordinate 
measured  eastward  from  tl»e  initial  longitude,  and  (I  is  the  true  latitude.  Here, 
as  in  the  wind  problem,  the  solution  is  resolved  into  two  independent  parts: 
tiu'  computation  of  (a)  a  mean  displacement,  /V,  and  (b)  a  deviation  displace- 
menl  with  components  /  ,((/'-  40).  Alsu,  A  =  J  v  +  /  .  As  was  pre¬ 
viously  stated,  any  displacement  in  latitude  involves  a  cliange  in  both  U  and 
'j.  riuis,  unlike  the  wind  problem  where  U  Is  uidepeiidenl  of  v  ,  in  this  ease, 
/\  is  a  function  of  ifl  . 
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For  reasons  which  v/ill  become  apparent  later,  it  is  necessary  to  solve 
the  deviation  part  of  the  problem  first.  Figure  4  is  a  solution,  for  40^N  lati¬ 
tude,  of  the  transformation  equations  (Eqs.  [4J  and  [S])  developed  in  Appen¬ 
dix  IV.  The  radiating  curves  are  arcs  of  great  circles,  m.arked  off  in  values 
of  s,  the  linear  distance  from  the  origin.  Each  great  circle  is  identified  by  a, 

I 

a  spherical  angle  measured  from  an  arbitrary  meridian.  Note  that  both  /  and 
0  depend  a  and  s.  Thus  we  may  write  the  following  symbolic  relations. 
(Note:  The  actual  relations  are  Eqs.  [4]  and  [S]  of  Appendix  IV.) 

A  =  A  (a,  s) 

0  -  B  (a,  s) 

The  general  method  of  solving  the  deviation  part  of  the  problem  consists  of: 

(a)  Integrating  Eq,  (11)  with  respect  to  s  along  each  value  of  a. 

(b)  Determining  S^,  the  appropriate  value  of  s  for  P  e  0.  50, 

(c)  Finding  the  A  and  0  values  corresponding  to  S  ,  using  Fig¬ 
ure  4, or  Eqs.  (4)  and  (5)  of  Appendix  IV. 

Before  proceeding,  however,  we  must  first  decide  on  a  metliod  of  accounting 
for  the  variation  of  6  with  latitude.  We  shall,  therefore,  describe  two  dif¬ 
ferent  methods  of  solution  and  compare  the  results  obtained  from  each  for  this 
particular  problem.  The  first  metliod,  called  tlie  "simplified  solution,  "  is 
quasi-linuar,  while  the  second  method  or  "complete  solution"  is  non-linear. 
From  a  computational  standpoint,  the  essential  difference  between  tlie  two  meth¬ 
ods  is;  In  the  simplified  solution  the  d  profile  is  introduced  after  integration 
of  tile  probability  function,  wliereas  in  the  complete  solution  the  d  profile  enters 
the  problem  before  integration. 

We  now  introduce  the  circumflex  symbol  (  )  to  denote  a  quantity  wiiich 

is  averaged  between  the  initial  latitude  (in  this  case  40^\\)  and  latitiuh;  0  . 

It  will  be  uiulersiootl  that  the  melhud  of  averaging  is  to  be  specified  in  eaclicase. 


l,i 


LONGITUDE) 


and 


S  =  d  (in  4)^/^  =  G  (1)  (3.  (in  4)^/^  . 

and,  since  G  (1)  =  0.86', 

S  =  0.  86  X  1.  176  d.  =  1.01  d.  (13) 

where  0\  is  expressed  in  degi'eea  latitude  per  day.  Since  S  is  independent 
of  <1,  the  solution  of  the  deviation  displacements  for  constant  mean  wind  is  a 
circle  on  the  earth's  surface  with  a  radius  which  is  a  function  of  U'.  If  we 
compare  Eqs.  (13)  and  (13a),  also  Eq.  (14a),  we  see  that  S  is  directly  pro¬ 
portional  to  C .  In  fact,  S  =  0.  86C,  We  have  seen  that  C  represents  the  ra- 
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dius  of  the  wind  deviation  circle.  However,  the  radius  of  the  displacement 
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deviation  jircle  is  not  S,  but  ■  sin  S.  This  difference  aris'^s  because  the 
radius  lies  in  the  plane  of  the  circle  and  not  on  the  curved  surface  of  the  earth 
where  S  lies.  Figure  5  is  an  enlargement  of  Figure  4  on  which  the  deviation 
displacements  of  .the  0.50  probability  curve  have  been  plotted  for  various  values 
of  U,  Since  these  curves  are  axially  symmetrical,  only  half  of  the  figure  is 
shown. 

At  this  point  v/o  introduce  the  profile  of  U  .  Wo  assume  that  the  solution 
at  latitude,  0  ,  is  the  same  as  the  solution  for  a  constant  mean  wind  of  mag> 
nitudu,  U  ,  whore 


The  protilc  of  U  is  shown  in  Figure  2,  page  11,  The  mapping  of  the  0,50  prob¬ 
ability  curve  for  the  deviation  part  of  the  solution  is  shown  in  Figure  5. 


0.  COMPI  ETE  SOLUTION  OF  THE  LISFLA CEMENT  DEVIATIONS 
Equation  (II)  is  rewritten: 


p  («)  =  exp 
d^- 


1 

id"  ! 


(lib) 


15 


x' (DEVIATION  LONGITUDE) 


“i ,  SnUj'.  H'li’.S  "rtf  ihr  <S 


di spS;\v  f'uu'rtl  0.^0  pJMh.ihUriV 


corvi' 


,m 


d  is  a  value  of  d  appropriate  for  the  population  of  displacement  deviations 
from  0  to  S.  Within  the  integration,  d  is  a  constant,*  Integrating  Eq.  (lib); 

,.S 


P  (a.S) 


i 


/ 

% 


8  exp 


26^  i 


ds 


=  I  1  -  exp 


16‘ 


(12b) 


a 


6  (  in  4)^^^  I  =  [ G  (1)0.  (  in  4)^'^ 


(13b) 


ct 

r.S 


i.Ol  0. 


Ja 


W  e  define  8.  =  J’  ds/^  . 


(14b) 


'm 


In  general,  0^  will  vary  with  a.  Thua,  unlike  the  solution  for  a  conatant 
mean  wind,  the  complete  solution  is  not  circular. 


*  The  reason  for  thia  atema  from  the  definition  of  a  probability  distribution 
i  -  ction.  To  show  thisi  we  will  integrate  Eq.  (U)  in  a  straightforward  man¬ 
ner,  considering  0  =  0(8): 


P(S)  s  / 


exp 


Integrating  by  parts: 

P(S)  =  1  -  exp 


s', 


‘20 


2/ 


da  . 


.03 


[  2fi"  1  6^ 


exp 


A  a  S  I 


,•0 


P(S)  -►  1  + 


-^0, 


S  ^2 
0-' 


exp 


20^  ^ 


20 


<2 


I  d0 


'  d0 


«i] 


iv'-;. 


!■ 


By  definition,  P  =  1  at  S  =  *  ;  but,  in  general,  the  integral  term  on  the  right 
does  not  vanish.  Therefore,  the  above  equation  ia  not  a  probability  function. 
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Table  3  shows  the  results  for  this  example,  obtained  by  numerical 
integration. 


Table  3.  Computation  of  displacement  deviation  (complete  solution) 


;  DEGREES 
LONGITUDE, 
ROTATED 
SYSTEM 

6. 

1 

DEGREES 
LA  TITU  DE 

DEGREES 

LONGITUDE 

a 

S 

0 

i 

/ 

0 

13.  33 

13.  46 

26.  5 

.  32.5 

45 

14.  08 

14.  22 

29.4 

36.  2 

90 

14.  61 

14.  76 

38.  5 

40.  2 

135 

12.40 

12.  52 

47.4 

32.  2 

180 

11.  65 

11.  87 

51.  9 

30.  0 

225 

12.39 

12.52 

47.4 

-32.  2 

270 

24.  61 

14.  76 

38.  5 

-40.  2 

315 

14.  08 

14.22 

29.4 

-36.  2 

The  solution  is  plotted  in  Figure  5.  The  difference  in  results  from  Uie  two 
methods  for  the  deviation  part  of  the  solution  is  practically  negligible. 


7.  SOLUTION  OF  THF  Ml’,h\N  DISPLACEMENT 

The  mean  longitude  displacement  for  a  particle  moving  along  a  latitude 
circle  with  the  speed  of  the  n\ean  wind  fur  that  latitude  is 

/t=  (U  sec  D)  T  . 

In  the  cvjmplete  solution,  where  we  consider  displacements  along  groat  circles 


\  -  (U  sec  1^)^  T 


(15b) 


w  h  e  r 


/  a 


/tt  .1,..*  n  \ 


0  (U  sec  d  )  ds  1 


In  the  simplified,  solution 


7\=  U  (sec  0')  T  (16b) 

a 

r0  ^  J'^a 

where  again  U  =  - - and  (sec  0)  =  °  -- 

0-40  S 

a 

The  two  solutions  are  plotted  in  Figure  6.  It  now  remains  to  add  the  two 
parts  of  the  solution.  The  final  results  are  shown  in  Figure  7.  It  may  be  seen 
that  the  difference  between  the  two  solutions  is  quite  small.  Nearly  all  of  the 
difference  is  produced  by  errors  in  the  simplified  solution  of  the  mean  displace¬ 
ment,  Since  the  mean  displacement  is  proportional  to  tirne,  there  wifi  be  a 
tendency  for  the  errors  to  increase  with  time.  On  the  other  hand,  the  fact  that 
S  also  increases  with  time  will  tend  to  diminish  the  errors.  On  the  whole,  it 

Q 

may  be  said  that  the  simplified  solution  yields  a  good  approximation. 

8.  MAPPING  OF  THE  PROBABILITY  FIELD 

To  illustrate  the  nature  of  the  probability  field,  we  have  plotted  probabil¬ 
ity  maps  for  an  initial  point  of  40°N  and  120°W  at  50,  000  feet  for  January  and 
July,  and  for  time  periods  of  1,  2,  and  3  days  (Figures  8  through  13).  The 
computation  method  used  to  obtain  this  series  of  charts  was  precisely  tlie  same 
as  outlined  for  the  cojnplete  solution.  These  maps  may  be  interpreted  or  used 
in  several  ways.  As  a  simple  example,  suppose  tl\at  a  constant-level  balloon 
is  lavinchod  at  50,000  feet  in  January  from  the  initial  point.  If  we  have  no  fur¬ 
ther  information,  then  the  indicated  probability  gives  the  correct  odds  that  tl\c 
balloon  will  be  located  williin  the  area  enclosed  by  the  given  probability  curve 
after  t  days.  Thus,  the  probability  is  0.25  that  the  balloon  is  within  the  0.  25 
curve  versus  0.75  that  it  i:  outside.  Similarly,  llu;  probability  is  0,99  tl\at 
the  balloon  is  witliin  the  0.99  curve  versus  0.01  that  it  lies  outside.  ‘I'he  moat 
probable  area  is  enclosed  by  the  0.50  probability  curve. 

With  suilaViltf  refinements,  this  presentation  could  possibly  he  used  for 
problems  involving  atomic  fallout.  In  tins  case,  we  would  need  dislrilnilion 
lunciivins  at  several  almosplieric.  levels.  We  cannot,  of  course,  proviile  detailed 
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information  on  the  amount  of  concentration  since  this  is  a  small-scale  phen¬ 
omena,  However,  we  can  provide  information  on  the  location  of  high  concen¬ 
tration  areas. 

The  foregoing  probability  statements  were  made  with  the  understanding 
that  our  only  available  information  is  the  time  and  place  of  the  balloon  launch¬ 
ing.  Now,  suppose  that  we  are  given  one  further  bit  of  information,  namely, 
the  initial  wind  velocity.  Our  theory  permits  us  to  construct  a  different  set  of 
probability  maps  which  take  into  account  this  information.  In  the  first  few  days, 
the  areas  enclosing  the  probability  curves  would  be  smaller  and  the  central  path 
would  be  different.  With  the  passage  of  time,  differences  between  the  two  sets 
of  maps  would  tend  to  disappear.  In  other  words,  the  added  information  re¬ 
duces  some  of  the  uncertainty  about  the  balloon's  location,  but  the  information 
becomes  less  important  with  increasing  time. 


9.  THE  DOWNSTKflAM  PROBABILITY  FUNCTION 


For  balloon  operations,  the  most  useful  probability  statement. is  the 
"Downstream  Probability  Density  Function"  (DPDF).  This  may  be  d’efined  as 
the  probability  density  with  respect  to  latitude  of  an  x-clisplacement  at  least 
as  large  as  a  specified  value  and  occurring  within  T  days.  Mathematically, 
this  density  function  is  expressed  as  follows: 


DPDF 


If  we  wish  to  interpret  tins  function  in  terms  of  balloon  operations,  we 

first  assume  that  a  large  lum-iluT  (,)f  balloons  are  launcliud  from  a  single  j>uinl 

over  a  long  period  of  time.  In  tins  ease,  x  and  y  are  eompont:nts  measured 

froiti  the  launching  site  anti  l  refers  to  the  time  interval  after  launching.  Thu 

term  (I  -  P  )  represents  the  (iroportional  number  of  balloons  to  be  found 

downstreaivt  from  x  at  lime  t  ;  '  •  (1  -  I*  )  rei)resents  tlu;  rale  at  which  bal- 

1  ^  ‘ 

loons  (low  past  x  at  lime  t,  p  :  (1  -  P_  )  represents  the  rale  at  which 

y » I  t  »  V 

b.Uloon.s  flow  downstream  past  tlie  point  x,  y  at  time  t.  l-'mally,  Kq.  (IV) 
represents  llu-  proportional  number  i>f  balloons  wlucli  may  be  exptu'ted  to  have 


.18 


flowed  downstream  past  the  point  x,  y  from  the  time  of  launching  to  the  time  T  „ 

The  various  steps  in  the  computations  required  to  solve  Eq,  (17)  for  a 

specific  example  are  illustrated  in  Figure  14,  For  convenience,  we  substitute 

for  p  the  j:vobabiiity  for  a  small  increment  of  y,  namely,  Ap  t  ‘ 

y ,  t.  '  y  >  ^  7  ‘  y  y 

The  DPDF  musi  be  mapped  on  a  spherical  earth  and  the  variation  of  6 
and  U  with  latitude  must  be  taken  into  account,  just  as  in  the  case  of  the  polar 
integration.  It  has  been  previously  demonstrated  that  use  of  the  simplified  so¬ 
lution  introduces  very  little  error  into  the  results.  For  this  reason  it  was  de¬ 
cided  to  adopt  the  simplified  solution  for  computing  the  DPDF.  Accordingly, 

Eq.  (17)  was  solved  for  various  values  of  U,  s,  y,  and  t.  The  next  step  is  to 
eliminate  x  by  the  transformation,  x  sec  p.  For  y  we  substitute. Q- 
We  now  have  a  table  with  the  parameters,  U,  X  ,  and  P.  This  is  Table  I, 
Volume  li.  The  final  step  is  to  substitute  the  values  of  U  for  U  for  each  re¬ 
gion,  season,  and  height.  Here  again 


U  = 


U  d  0 


VV  i:  iio\s'  have  a  set  of  tables  for  each  region,  season,  and  luhght  tor  various 
value's  of  (1  and  p.  Two  such  tables  appear  in  Volume  III.  namely,  Table  1, 

-North  '  morica,  and  Table  II,  rdurasia. 

1  lu  re  iro  two  limitations  on  the  downstroatn  probability  concept.  In 
the  fir.st  place,  on  a  spl.e-rical  earth  some  air  particles  may  take  a  decided 
'lorthward  cotinie  and  reappear  at  a  lower  latitude  very  far  downstream.  While 
■.V ..tii d  In;  feasible  to  formulaic  a  solution  wliich  incorporates  this  type  of  tra- 
i." u  wi;  decided  to  exclude  these  cases  from  consideration.  Tbu-s  the  com- 
pu'.  ;,  ;  d  ir.i  ii  Volumes  II  and  III  refer  only  to  liiose  motions  which  are  cssen- 
‘  in!  U  (!>  ■■...  ■  I  m. 

!  he  Hecunu  1  i mita! ion  vu.curs  when  the  mean  vtonal  wind  ciiangcs  sign  at 
smjn  l.ttiPuio,  t'lus  occurs,  for  t-xample,  over  North  America  at  uO.tUIO  fe«-t 
•imeiner  and  ai  HO,  t't.U)  fe.  i  in  the  fall  for  certain  initial  l.ititudes.  Our  solu- 
h.i'i'e  IS  Eo  make  two  sets  of  tables,  oiu*  for  eastward  displacements  and 
d'.c  fdir  westward  viisplacenu-c.is. 
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To  obtain  a  coniplete  solution  of  the  wind  climatology  problem,  certain 
basic  quantities  must  be  plugged  imo  the  statistical  equations.  These  include 
the  varietion  by  season,  latitude,  and  height  of  liie  mean  instantaneous  wind 
and  its  associated  distribution  functions.  This  is  analogous  to  the  solution  to 
a  vast  jigsaw  puzzle  where  many  important  pieces  are  missing. 

The  first  ••■tep  in  rhe  solution  was  a  rather  complete  survey  of  available 
literature.  From  this  we  were  able  to  deduce  a  somewhat  sketchy  broad-scale 
picture  of  the  wind  distribution  at  upper  levels  over  the  earth.  It  became  im- 
meaialely  apuareni  that  fur  more  details  \v  <■:.  vvoulcl  hs-vc  to  resort  to  indirect 
or  "bootstrap”  techniques. 

One  of  the  (hasired  quantities  is  some  measure  of  the  distribution  asso¬ 
ciated  with  the  m.ean  wind.  If  the  wind  is  normally  distributed,  a  single  para¬ 
meter,  such  as  the  standard  vector  deviation  or  the  constancy  (ratio  between 
the  vector  mean  wind  and  the  .scalar  mean  wind),  can  be  used.  ^  For  conven¬ 
ience  and  ease  of  computation,  it  was  decided  to  compute  the  wind  constancy. 
To  minimize  the  known  bias  in  the  observed  upper  wind  records,  the  differ¬ 
ence  method  of  averaging  w.ss  used  throughout  the  study. 

We  decided  to  concentruto  on  an  intensive  analysis  of  ihnu'  data  sou.rcQs; 
(a)  a  9-year  record  of  15  U.  S.  Weather  Bureau  Stations;  (b)  a  .senes  of  ^  to 
4  years  ul  upper  air  synoptic  uiap's  of  the  West  cicrnuin  .Vletcorologicttl  Ser- 

g 

vice  at  i  levels:  500  mb  (20,  i)!)0  ft).  225  mb  (36,000  ft)  and  96  mb  (54,  000  ft): 
and  (c)  rawuisonde  records  from  Berlin  and  .several  other  Gernutn  sialions. 

kach  ot  llie  German  syavii'lic  chai'ts  for  the  months  of  .tanuarv.  April, 
July  and  October  was  pliolographical ly  enlarged  to  a  size  of  15”  by  20”,  and 
isotachs  (lines  ot  tu'iual  wiml  speed)  were  constructed,  usin;.*  the  computed 
geostrophic  winds  as  welt  as  obs er ved  \.v  ind  speeds.  l''roni  tiu-;  analyzed  eharts, 
mean  vector  and  scalar  winds  svere  eonqniled  for  eai’li  latitiule,  season,  and 
height,  combining  the  data  frop'  i.nigitude  30'\v  to  JO^’k  into  a  .single  statistic. 


In  all,  1000  maps  were  3,nalyzed  and  35  wind  valuns  coniputed  from  each  map. 

The  United  States  data  were  summarized  by  punciied  cards.  The  ratio 
between  the  scalar  winds  a.t  any  level  ai\d  some  standard  level  v/as  the  most 
stable  statistic  that  could  be  found.  For  this  reason  a!'  vi'iC  da.ta  were  re¬ 
duced  to  ratios  between  the  scalar  wind  speeds  at  given  levels  and  the  cor- 
respoTiding  scalar  winds  ar  b  :ci.lometers,  oinoniarizcd  values  of  this  sta¬ 
tistic  were  smoothed  by  hirmonic  analysis  ae(i  matched.  for  level,  with 

similar  data  from.  EuropeS’i  sources.  In  this  way  the  latitudinal  variation  of 
the  scalar  ratio  vva.s  established.  For  levels  above  9b  n'.b  ib,.’  Gorman  wind 
data  were  used  in  the  comparison.  These  were  supplemented  by  data  from 
nuannscript  Nori'nern  i-i^-mispher e  sytioptic  maps  lor  1953  at  100  mb  (53,000  ft) 
and  50  mb  (68,000  tt).  higure  15  is  an  exanipie  of  one  .sol.  of  ratio  analysis. 

The  trough  lines  in  the  upper  portion  of  the  chart  delineate  the  boundaries  of 

<  Vl  fc.’i  <  •  f  l.'N  1  1  J. .  U 

Maving  established  tlm  scalar  wiiid  ratios,  it  was  tlien  possible  to  com¬ 
pute  reliable  scalar  wind  speed.s  by  muliiplying  the  ’’atios  with  the  appropriate 
mean  .scalar  winds  al  '^00  mb  which  are  known  wiiii  ,,i;- 1 .  ic  bU;  accuracy. 

This  left  the  vector  mt-an  winds  to  he  dei  "'-m  [t;  was  found  that  this 

stalislic  was  diificult  to  detcrimae  ac.curatel''  :  -  'P'  a  limited  sample  because 
of  variations  in  wind  direction.  An  intensive  .studv  of  the  500-mb  United  States 
data  led  to  the  discoi  erv  that  a  remarkably  men i;--!  rut  relationship  existed  be¬ 
tween  the  scalar  and  \'ee'>'i’  > ici n  inds .  It  u-  'a  e* ;  known  flial  strong  winds 
are  more  eonstant  than  wmil  vinds.  thn  ’.uth'  Lietween  the  vector  and 

actilar  'oni  i^  a  me.isure  te  -  n.d  coi'staucy,  t’ais  ratio  should  rise  as  the  wind 
me  r  cM ,s es .  I'his  a -is  :\>',ni(i  in  t)e  tlie  ease.  It  was  ,il.so  foinut  tlial  the  presence 
of  the  mean  Ji't  .'hreum  pre>duees  an  ■uhiitiun.'U  mcreatie  in  constancy  over  and 
aluivi  :ie  eltect  of  wind  si->eeil  .ilone,  1  hm  eii.'ct  m  also  known  qu.il  italively, 
‘i'hese  rel.it  ionships  .in'  shown  in  higure  Je,  1  In-  degree  t.e  w-hu'h  'he  United 
.St.ites  ''I'lt-mh  rl.ita  tUtid  tU,  ife  ri  veil  car',  i-'  'C.  .  id  the  p''s  s  ibilit  v  that  this 

was  a  universal  .statistu  .il  l.i'.v  governing  the  w  e  ,a  i  c- r  1 1  es  .  A  1 1 1'tupl  ing  to  verify 
thi.s  su g esiu'c,  siriagenil;.,  ,ve  u.sed  test  da'.i  from  .■  d.iai.mi  regioji,  ,a  rniuoto 
level,  and  lor  a  ui’icreni  time  '.leruui.  n.im  -lv,  th--  Mi'i.i  \',‘md  d.ita  fer  Europe, 
lahli’  •!  .‘i'n.'ws  th.i’  esv  .•llent  rx'.sulis  o!  this  te.st,  > 'vusimuentlv,  we  may  conclude 


that  the  mean  wind  uniquely  determines  its  distribution  function. 


Table  4,  Correlation  of  teat  data  (96-mb  winds  kts  .  over  Europe) 


MONTH 

LATITUDE 

SCALAR  WIND 

VECTOR 

Predicted 

WIND 

Observed 

Jan 

70 

37.  0 

27.  2 

26.9 

Jan 

60 

46.  3 

37.4 

36.9 

Jan 

50 

37.4 

27.  5 

26.9 

Jan 

40 

37.  6 

27.  6 

27.  0 

Jan 

35 

40.9 

32.  I 

31.  3 

Apr 

70 

27.  3 

19.  3 

17.  1 

Apr 

60 

28.  8 

20.6 

21.4 

Apr 

50 

22.  9 

15.  3 

13.4 

Apr 

40 

29.  0 

20.7  . 

19.9 

Apr 

35 

35.  7 

27.  0 

•  28.3 

Jul 

70 

16.  9 

8.5 

7.0 

Jul 

60 

20.  0 

12.4 

11.6 

Jul 

50 

20,  6 

13,  0 

11.0 

Jul 

40 

24.  8 

17.2 

r6.7  1 

Jul 

35 

26.  8 

18.9 

21.  0 

Oct 

70 

27.  5 

19.5 

i9.  8 

Oct 

•  60 

28.  6 

20.4 

19.8 

Oct 

50  ■ 

25,  3 

17.  6 

16.0 

Oct 

40 

27.  5 

19.5 

18.0  1 

Oct 

35 

29.7 

21.2 

21.9  ' 

I 


Correlation  coefficient,  r  =0.99  i 

A  similar  relationship  was  derived  for  the  easterlies  from  what  little  data 
were  available,  Considerably  less  confidence  can  be  placed  in  these  results. 
Nevertheless,  indications  are  that  the  high-level  easterlies  are  extremely  con¬ 
stant  winds. 

Using  the  data  from  Figure  16  we  computed  the  relation  between  the  vec¬ 
tor  moan  wind  and  the  standard  deviation  (Figure  17), 


APPENDIX  II 


UPPER  WINDS  OVER  EUROPE  AND  NORTH  AMERICA 


The  mean  zonal  winds  and  standard  deviations  used  in  this  study  are 
presented  in  Table  5.  Above  54,000  leet  the  sparse  European  data. was  sup¬ 
plemented  by  more  complete  data  from  North  Am.erica„  Thus  the  winds  at 
high  levels  in  Europe,  while  not  as  reliable  as  the  lower  winds,  m.ay  be  re¬ 
garded  as  a  first  approximation  to  the  true  winds.  The  standard  deviations 
in  Table  5  are  taken  from  Figure  17,  Appendix  I. 

Figures  18  through  23  summarize  the  important  variations  of  the  winds 
in  Table  5,  namely,  seasonal,  latitudinal,  and  vertical.  The  80,  000 -foot 
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are  characteristic  of  the  layer  below  50,000  feet. 


1.  SEASONAL  VARIATION  (Figures  18  and  19} 

The  seasonal  variation  of  winds  at  40,  000  and  80,  000  feet  over  North 
Amei'ica  at  40  N  {Figure  18)  and  Evirope  at  50''^N'  (Figure  19)  was  obtained  by 
harmonic  analysis.  Over  botii  regions  a  siui’P'  wave  was  sufficient  to  repre¬ 
sent  the  80,  OOO-foot  data  while  two  waves  wore  rv<iuirod  to  fit  the  40,000-foot 
winds. 

lire  bO,  OOO-foot  seasonal  s’ariation  tor  Iwilh  continents  consists  of  a 
single  maximum  oi  v-vesterlies  in  midwinter  and  a  single  maximum  of  easter¬ 
lies  in  tnidsummer  with  reversals  of  direction  m  April  and  September. 

At  40,000  leet  tin''  doul)leAva\'e  strnctnre  iii  tlii  w I'stt.'t’l i es  is  very  pro¬ 
nounced  in  Europe  and  ijuite  small  in  Nortii  Ann'rica,  I'lu're  are  two  distinct 
maxinta  in  I'.urope,  the  primary  ui  Marc'n  aiui  the  seeondary  in  November;  and 
two  minima,  tlu'  majc'i*  in  August  .’.nd  nuiuir  in  I'el) nui ry .  lly'  contrast,  the 
North  A'inencan  curves  sliew  a  \ery  prunuunc ni  inaximum  of  100  kls  in  January 
and  a  mere  suggi'stion  oi  a  seiuiuLiry  maximum  in  Julv,  I'he  effect  of  tlte 
seeond  '.vave  in  lliis  c.'ase  is  to  tlalien  out  the  simnner  minimum  (lluit  is,  east¬ 
erly  maximum)  so  Unit  llu-  mean  r.en.ii  .i  ind  is  sirtnally  constant  from  .May 
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and  standard  deviation. 


through  September.  Note  that  the  North  Amerii'.cin  wmds  are  more  than  twice 
as  strong  as  the  European  winds  in  winter  at  llie  40,U00-foot  Level, 

2.  LATITUDINAL  VARIATION  (Figures  20  and  21) 

At  40,  000  feet  both  the  winter  and  summer  latitudinal  profiles  for  North 
America  each  show  a  single,  sharp  maximum  (Figure  20)  while  the  corres¬ 
ponding  curves  for  Europe  display  double  maxima  (Figure  21).  A  comparison 
of  the  January  and  July  curves  revc  ils  that  the  European  maxima  at  60*^N  and 
33°N  remain  fixed  in  position,  whereas  the  single  maximum  over  North  Amer¬ 
ica  migrates  northward  from  40°N  in  January  to  45°N  in  July,  Note  also  the 
difference  in  range  of  the  mean  zonal  w'inds  with  respect  to  latitude  i-n  North 
America  as  compared  with  Europe,  The  range  is  far  greater  in  North  Amer¬ 
ica  than  in  Europe,  being  more  than  three  times  as  great  in  January  and  in 
excess  of  twice  as  large  in  July. 

Of  considerable  interesL  is  tlie  transition  zone  in  the  vicinity  of  50,000 
feet  where  the  European  \vinds  change  from  a  double  maxima  pattern  below  to 
a  single  maximum  above,  while  the  North  American  profiles  change  in  pre¬ 
cisely  the  opposite  direction  from  a  single  ma.ximum  below  50,000  feet  to 
double  maxima  above. 

In  January  the  westerlies  at  80,  000  feet  over  Europe  reach  a  broad 

maximum  at  60^N,  while  in  Nortli  .America  the  major  maximum  is  north  of 

0  o 

70  N  and  there  is  a  secondary  peak  at  about  37  N.  In  .'uly  the  winds  are 

easterly  above  50,000  foot  at  aU  latitudes  over  both  iJurope  and  Nortli  Amer- 

o 

lea.  In  Europe  tlie  easterlies  vary  only  .stightlv  with  latitude  from  25  N  to 
75^\'.  There  is  a  flat  maximum  In  the  vicinUv  of  55'^N.  In  contrast,  the 
summer  higti-level  easterlies  over  Nortli  Anu  rica  have  two  maxima,  one 
north  of  75'’^N  and  t!u'  oUu'r  south  of  2-)^  As  m  the  case  of  the  winds  be¬ 
low  50,000  ft'ol.  tlu'  uppor-Uwel  a  uids  over  \oriU  Ameriivi  have  a  muchlarger 
range  than  the  European  Ainds,  tlie  r.uiee  at  tii’.uou  leet  being  niore  than  six 
times  as  great  in  winter  and  tt\e  tinu';,  nS  larg*  m  .'■uimmer. 


MEAN  ZONAL  WINDS  (KNOTS) 


F\,''ure  IB.  Seasoiial  variation  of  winds 
aver  North  America 
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Figure  20.  Latitudinal  variation  of  winds 
over  North  America 
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3.  VERTICAL  VARIATION  (Figures  22  and  23) 


The  variation  with  height  of  the  mean  zonal  wind  is  similar  for  North 
America  and  Europe.  The  mean  zonal  wind  reaches  a  maximum  at  40,  000 
feet  throughout  the  year  between  25°N  and  75°N.  In  January  the  wind  de¬ 
creases  with  height  above  40,  000  feet  to  a  minimum,  then  increases  again 
up  to  some  level,  presently  unknown.  Over  Europe  the  minimum  occurs 
at  70,  000  feet;  over  North  America,  at  80,  000  feet  or  above.  The  level  of 
minimum  wind  does  not  appear  to  vary  significantly  with  latitude.  In  July 
the  wind  also  .decreases  with  height  above  40,  000  feet.  In  the  vicinity  of 
60,  000  feet  the  di  rection  reverses  from  West  to  East.  Above  this  level 
the  easterlies  increase  with  height.  Over  North  America,  to  the  noBth  of 
70°N  and  to  the  south  of  30°N,  the  winds  are  easterly  at  all  levels  in  summer. 

In  middle  latitudes  in  January  the  westerlies  over  North  America  are 
stronger  at  all  levels  than  over  Europe.  The  greatest  difference  occurs  at 
40,000  feet,  the  level  of  maximum  wind,  where  a  speed  of  100  kts  occurs 
in  North  America  at  40°N  as  compared  with  45  kts  in  Europe  at  35°N. 

(Note- the  two  latitudes  selected  for  this  comparison  both  lie  along  or  near 
the  horizontal  axis  of  the  mean  wind  maximum  for  the  region  in  question. ) 

In  summer  it  appears  that  the  maximum  easterlies  occur  above  80,000 
feel  in  Nortli  America  and  near  80,  000  feet  in  Europe.  Of  special  interest  is 
the  fact  that  the  high-level  easterlies  have  the  greatest  constancy  of  all  the 
wind  regimes  we  have  studied. 


APPENDIX  III 


VARIATION  OF  THE  STANDARD  DEVIATION  WITH  TIME 


We  have  adopted  the  classical  diffusion  theory  developed  by  G.  I.  Tavlor. 

He  considered  a  fluid  containing  uniformly  distributed;  small-scale  eddies  which 
produce  diffusion  by  continuous  movement.  We  have  adapted  this  theory  to 
large-scale  eddies  of  the  general  circulation. 

If  R  is  the  correlation  between  u  (the  West  -  East  wind  component)  at 
time  t  and  u  at  time  t  +  f ,  where  ^  la  a  small  time  increment,  then  Taylor  shows 
that 


(x)  =  (u) 


R,  d^  dt 


(1) 


where  6  (x)  is  the  standard  deviation  of  x-diaplacementa  (that  is,  West  -  East 
displacements)  and  d  (u)  is  the  standard  deviation  of  the  wind. 

The  solution  of  this  equation  depends  npoji  the  form  of  the  auto-correlation 
function,  R^  .  Brunt^  derived  the  following  exponential  form  for  this  function  by 
considering  a  fluid  where  eddies  are  all  the  same  siae: 


H,  ^  e 


afc. 


where  a  is  a  constant  depending  on  the  sixe  of  the  eddies, 
Substituting  Eq.  (2)  into  Eq.  (1): 

•>  '•  f  ^ 


e  dr  dt  , 


a  a** 


U) 


(.3) 


or 


whore 


d(u) 


G  (T) 


G  (T) 


Vi) 


•5b 


We  have  seen  that  the  function  G(T),  which  defines  the  ratio  between  d(x) 

and  (5(u),  is  obtained  by  integrating  the  auto-correlation  function,  R  ,  twice 

% 

with  respect  to  time.  If  the  eddies  are  all  of  one  size  then  the  auto-correlation 
function  is  exponential  and  the  G-function  has  the  form  given  in  Eq,  '4a). 

Obviously,  eddies  in  nature  are  not  of  a  single  size;  therefore,  R  is  not 
strictly  exponential.  However,  the  principal  concern  here  is  not  the  exact  form 
of  R  ,  since  any  errors  in  this  function  are  smoothed  out  by  the  double  integra- 
tion  required  to  obtain  G(T),  The  important  consideration  is  the  degree  to  which 
Eq,  (4a)  approximates  the  actual  G-function. 

To  see  if  Eq.  (4a)  does  in  fact  closely  approximate  the  actual  G-function, 
we  computed  G(T)  directly  in  an  Eulerian  system  of  coordinates  (where  mea- 
surements  are  made  at  fixed  points  in  space  as  time  varies),.  The  advantage  of 
this  system  is  that  computations  are  made  from  actual  wind  data.  This  affords 
the  opportunity  of  determining  the  G-function  for  as  long  a  time  period  as  desired. 

Actually,  we  are  interested  in  the  G-function  for  a  Lagrangean  system 
(where  measurements  are  taken  along  the  path  of  a  moving  particle),  but  limited 
samples  of  trajectories  restrict  computation  of  the  function  to  periods  of  only 
2  or  3  days.  However,  there  is  no  reason  to  assume  that  the  form  of  the  G-func¬ 
tion  will  be  different  in  the  two  systems.  Therefore,  we  may  use  the  Eulerian 
computations  to  determine  the  form  of  the  G-function. 

The  Eulerian  computations  were  made  using  German  wind  data  at  50*^N 
for  225  millibars.  Winds  were  resolved  into  u  and  v  components  and  the  ratios 
and  were  computed  at  1-day  intervals  up  to  a  period  of  6  days.  It 

was  found  that  the  two  ratios  have  identical  time  variations,  thus  G^^(T)=Gy (T). 
Furthermore,  the  G-function  is  very  closely  approximated  by  Eq.  (4a)  with 
e  ^=0.6. 


While  there  is  every  reason  to  expect  that  the  form  of  the  G-function  is 
the  same  in  both  the  Eulerian  and  Lagrangean  systems,  it  should  not  bo  antici- 
pated  that  the  constant,  e  ,  is  the  same  in  the  two  systems.  In  fact,  we  would 
expect  this  constant  to  be  smaller  in  the  Lagi*angean  system,  since  space  as 
well  as  time  variations  are  included  in  the  ■'■omputation  of  G(T).  We  will  now 
discuss  tnree  separate  estimates  of  the  G-fuiiction  in  Lagrangean  coordinates. 
All  three  estimates  show  the  expected  smaller  value  of  the  constant,  e~^  . 
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A  sample  of  55  Moby  Dick  constant-level  balloon  trajectories  at  about 
200  mb  was  used  in  conjunction  with  associated  upper  air  charts  to  compute 
values  of  =  G(T)  for  time  intervals  of  12,  24,  and  36  hours.  A  value  of 
e  ^  =  0,4  was  obtained. 

Another  value  of  e  was  computed  from  trajectory  distributions  cal¬ 
culated  by  Spreen  (unpublished  report)  from  about  3000  geostrophic  trajecto¬ 
ries,  constructed  from  upper-level  charts  of  the  West  German  Weather  Ser- 
vice,.  Bad  Kissingen,  Germany.  The  value  e  =0.4  was  again  obtained. 

R.  R.  Rapp^^  computed  an  empirical  distribution  from  a  sample  of 

1 

179  trajectories  at  500  mb  for  two  fall  seasons.  Using  Taylor’s  development 
and  assuming  that  e  ^  =  0.5,  he  calculated  the  circular  normal  distribution 
for  the  same  3a,rnple.  The  X-  test  revealed  no  significant  difference  between 
the  empirical  and  theoretical  distributions;  therefore,  Rapp  concluded  that 
the  vdice  e  ^  =  0.5  ,  which  he  had  chosen,  was  reasonable. 

On  the  basis  of  the  above  studies  we  chose  e  ^  =  0.4  as  being  the  most 
representative  value.  Substituting  this  value  into  Eq.  (3)  we  obtain: 


6  (x.  T)  =  \/Jd  (u) 


1.  0917T  -  1.1919  (1  -  [0.4] 


(5) 


This  equation  specifies  the  variation  of  d(x)  with  time. 


APPENDIX  IV 


SPHERICAL  TRANSFORMATION 


Consider  a  coordinate  system  x,  y,  and  z  on  a  spherical  earth.  This 
coordinate  system  is  rotated  relative  to  the  conventional  coordinates  x,  y, 
and  z  in  such  a  way  that  the  x  and  z  axes  are  both  inclined  by  an  angle  0  = 
0^  to  the  X  and  z  axes,  respectively,  while  the  y  axis  coincides  with 
the  y  axis.  (See  Figure  24.)  The  angle,  0^,  is  a  fixed  latitude  on  the  real 
earth. 


Z  Z 


Figure  25.  Axes  s  and  z 


Any  arbitrary  poiat.  P  ,  oji  the  sphere  li.is  coordinates  y.  and  z  in 
the  conventional  system  and  x,  y,  and  z  in  the  rotated  system.  We  wish  to 
find  the  r elativ>nship  btuween  these  two  coordinate  systems,  Since  y  s  y  ,  we 
already  have  one  such  relaiuutslup,  hem  i-  the  problent  reduces  to  two  dimen¬ 
sions.  (See  Figure  ,1S.  )  Thus, 

X  =  X  sin  0  +  a  cos  0 

o  a 

y  =  y  ( 1 ) 

z  =  z  sin  0  -  X  cos  0 

o  o 


Inirodui’ing  the  fcdlowing  notation; 

K.  radius  of  the  ea;th. 

.Angles  (degr«*es);  Rotated  System 

s,  co-l.ililvuie 
‘t,  longitude 


C * -n V er\{ » una I  S V s t em 
0,  lalilude 
,  longitude 
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In  the  rotated  coordinate  system: 

X  =  R  sin  s  cos  a 
y  =  R  sin  s  sin  a 

z  =  R  cos  s  (Z) 

In  the  conventional  coordinate  system; 

X  =  R  cos  0  cos  A 
y  =  R  cos  0  sin  X 

z  =  R  sin  0  (3) 

Eliminate  x,  y,  z  and  x,  y,  z  from  Eq,  (1)  by  substituting  Eqs.  (Z) 
and  (3)  to  obtain  0  and  X  in  terms  of  0^,  a,  and  a. 

The  final  result  is; 

cos^  0  =  ain^  s  1  -  cos^  0^  (1  +  coa^  a)  | 

+  cos  0  +1^  sin  Zs  cos  Z  0  cos  (4) 

o  ^  o  a 

ctn  A  =  ctn  a  sin  0  +  ctn  s  cos  0  cac  a 
o  o 


■(5) 
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No.  14,  Atmospheric  r’Ku\  IMtierns  and  Their  Representation  by  SphericaUSurface  Harmonics,  II.  Huur* 
witz  and  Richard  A.  Craig,  Geophysics  Research  Directorate,  July  1952. 

No,  15.  HatrL-''catiering  of  Electromagnetic  Waves  From  Spheres  nntl  Spherical  Shells,  A,  I,.  Aden, 
tieophvsics  Research  Directorate,  July  1952. 

No.  In.  Notes  on  the  Theory  of  la\rge-Scale  Distnrbnnces  in  AiniiKpherie  Flow  With  Applications  to 
Nnmerii  al  Weather  Prediction,  I’hilip  Duncan  Thompson,  Major,  F.  S.  Air  l‘We,  (•eophysics 
Researcli  Directorate,  July  19.52. 
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GEOPHYSICAL  RESEARCH  PAPERS  (Continued) 


No.  17.  The  Observed  Mean  Field  of  Motion  of  the  Atmosphere,  Yale  Mintz  and  Gordon  Dean,  Geophysics 
Research  Directorate,  August  1952. 

No.  18.  The  Distribution  of  Radiational  Temperature  Change  in  the  Northern  Hemisphere  During  March, 
Julius  London,  Geophysics  Research  Directorate,  December  1952. 

No.  19.  International  Symposium  on  Atmospheric  Turbulence  in  the  Boundary  Layer,  Massachusetts  Insti¬ 
tute  of  Technology,  4-8  June  1951,  edited  by  E.  W.  Hewson,  Geophysics  Research  Directorate, 
December  1952. 

No.  20.  On  the  Phenomenon  of  the  Colored  Sun,  Especially  the  “Blue”  Sun  of  September  1950,  Rudolf 
Penndorf,  Geophysics  Research  Directorate,  April  1953. 

No.  21.  Absorption  Coefficients  of  Several  Atmospheric  Gases,  K..  Watanabe,  Murray  Zelikoff  and  Edward 
C,  Y.  Inn,  Geo,  -sics  Research  Directorate,  June  1953. 

No.  22.  Asymptotic  Approximation  for  the  Elastic  Normal  Modes  in  a  Stratified  Solid  Medium,  Norman  A. 
Haskell,  Geophysics  Research  Directorate,  August  1953. 

No.  23.  Forecasting  Relationships  Between  Upper  Level  Flow  and  Surface  Meteorological  Processes, 
J.  J.  George,  R.  0.  Roche,  H.  B.  Visacher,  R.  J.  Shafer,  P,  W.  Funke,  W.  R.  Diggers  and  R.  M, 
Whiting,  Geophysics  Research  Directorate,  August  1953. 

No.  24.  Contributions  to  the  Study  of  Planetary  Atmospheric  Circulations,  edited  by  Robert  M.  White, 
Geophysics  Research  Directorate,  November  1953. 

No.  25.  The  Vertical  Distribution  of  Mie  Particles  in  the  Troposphere,  R.  Penndorf,  Geophysics  Re¬ 
search  Directorate,  Mtrch  1954. 

No.  26.  Study  of  Atmospheric  Ions  in  a  Nonequilibrium  System,  C.  G.  Storgis,  Geophysics  Research 
Direclorute,  April  1954. 

No.  27.  Investigation  of  Micfobarometric  Oscillations  in  Eastern  Massachusetts,  E.  A,  ^'lau^nud,  A.  U. 
Mears,  F.  A.  Crowley,  Jr.,  and  A.  P.  Crary,  Geophysics  Rosenreh  Directorate,  \h»y  1954. 

No.  28.  The  Rotulion-Vibrntion  Spectra  of  Ammonia  in  the  6-  and  10-Micron  Regions,  R.  G.  Dreene,  Jr., 
Capt,,  U.SAF,  Geophysics  Research  Directwate,  June  1954. 

No.  29.  Seasonal  Trends  of  remperalure,  Density,  and  Prt^ssure  in  the  Stratosphere  Obtained  With  the 
SearcliUghi  Probing  Technique,  Louis  ElteruMn,  July  1954. 

No.  30.  Proceedings  of  the  Conference  on  Auroral  iMtysics,  edited  by  N.  C.  Gerson,  Geophysics  Re¬ 
search  Direclorute,  July  1954. 

No.  31.  Fog  Modification  by  tiold-Watcr  Seeding,  Vernon  (*.  Plnuh,  Geophysics  Research  DirecliMiile, 
August  1954. 


GEOPHYSICAL  RESEARCH  PAPERS  (Continued) 


No.  32.  Adsorption  Studies  of  Heterogeneous  Phase  Transitions,  S.  J.  Birstein,  Geophysics  Research 
Directorate,  December  1954. 

No.  33.  The  Latitudinal  and  Seasonal  Variations  of  the  Absorption  of  Solar  Radiation  by  Ozone, 
J.  Pressman,  Geophysics  Research  Directorate,  December  1954. 

No.  34.  Synoptic  Analysis  of  Convection  in  a  Rotating  Cylinder,  D.  Fultz  and  J.  Corn,  Geophysics 
Research  Directorate,  January  1955. 

No.  35.  Balance  Requirements  of  the  General  Circulation,  V.  P.  Starr  and  R.  M.  White,  Geophysics 
Research  Directorate,  December  1954. 

No.  36.  The  Mean  Molecular  Weight  of  the  Upper  Atmosphere,  Warren  E.  Thompson,  Geophysics  Re¬ 
search  Directorate,  May  1955. 

No.  37.  Proceedings  on  the  Conference  on  Interfacial  Phenomena  and  Nucleation. 

I.  Conference  on  Nucleation. 

II.  Conference  on  Nucleation  and  Surface  Tension. 

III.  Conference  on  Adsorption. 

Edited  by  H.  Reiss,  Geophysics  Research  Directorate,  July  1955. 

No.  38.  The  Stability  of  a  Simple  Baroclinic  Flow  With  Horizontal  Shear,  Leon  S.  Pocinki,  Geophysics 
Research  Directorate,  July  1955. 

No.  39.  The  Chemistry  and  Vertical  Distribution  of  the  Oxides  of  Nitrogen  in  the  Atmosphere,  L. 
Miller,  Geophysics  Research  Directorate,  April  1955. 

No.  40.  Near  Infrared  Transmission  Through  Synthetic  Atmospheres,  J.  N.  Howard,  Geophysics  Res- 
search  Directorate,  November  1955. 

No.  41.  The  Shift  and  Shape  of  Spectral  Lines,  R.  G.  Dreeno,  Geophysics  Research  Directorate, 
October  1955. 

No.  42.  Proceedings  on  the  (ionference  on  Atmospheric  Electricity,  R.  Holzer,  W.  Smith,  Geojdtysics 
Research  Directorate,  December  1955. 

No.  43.  Methwls  and  Results  of  Upper  Atmospheric  Research,  J.  Kaplan,  G.  Schilling,  H.  Kallnmn, 
Geophysics  Research  Directorate,  November  1955. 

No.  44.  Luminous  nnd  Spectral  Reflectance  us  Well  as  Golors  of  Natural  Objects,  R.  Penndorf,  Geo* 
physics  Research  Direclwate,  February  1056. 

No.  45.  New  Tables  of  Mie  Scattering  Functions  for  Spherical  Particles,  R.  Penndorf,  B.  Goldberg, 
Geophysics  Research  Dir«cl«>rate,  March  195ti. 

Nu.  46.  Results  of  Numerical  Forecasting  With  the  Uarolropic  and  'Uhermutropic  Mmlcls,  tt.  Gales, 
L.  S.  Pocinki.  b.  Jenkins,  Geophysics  Research  Directorate,  April  1956. 
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No,  47.  A  Meteorological  Analysis  of  Clear  Air  Turbulence  (A  Report  on  the  U.  S.  Synoptic  High- 
Altitude  Gust  Program),  H.  Lake,  Geophysics  Research  Directorate,  February  1956. 

No.  48.  A  Review  of  Charge  Transfer  Processes  in  Gases,  S.  N.  Ghosh,  W.  F.  Sheridan,  J.  A.  Dillon, 
Jr.,  and  IL  D.  Edwards,  Geophysics  Research  Directorate,  July  1955. 

No.  49.  Theory  of  Motion  of  a  Thin  Metallic  Cylinder  Carrying  a  High  Current,  C.  W.  Dubs,  Geo¬ 
physics  Research  Directorate,  October  1955. 

Hurricane  Edna,  1954:  Analysis  of  Radar,  Aircraft,  and  Synoptic  Data,  E.  Kessler,  III  and 
1).  Atlas,  Geophysics  Research  Directorate,  July  1956. 

Cloud  Refractive  Index  Studies,  R.  M.  Cunningham,  V.  G.  Plank,  and  C.  F.  Campen,  Jr. 
Geophysics  Research  Directorate,  October  1956. 

A  Meteorological  Study  of  Radar  Angels,  V.  G.  Plank,  Geophysics  Research  Directorate, 
Viiutist  1956. 

No.  Tvl.  riu'  Construction  and  I’se  of  Forecast  Rejisters,  1.  Gringorten,  1.  Lund,  M.  Miller,  Geo¬ 
physics  Research  Directorate,  June,  1956. 

No.  51.  Solar  (.eoniagnetic  Ionospheric  Parameters  as  Indices  of  Solar  Activity,  F.  Ward  Jr.,  Goo- 
plivsics  Research  Directorate,  November,  1956. 

No,  55.  Preparation  of  Mutually  Consistent  Magnetic  Charts,  Paul  Fougoro,  J.  McClay,  Geophysics 
Research  Directorate,  June  1957. 

No,  56.  Radar  Synoptic  Analysis  of  an  Intense  Winter  Storm,  Edwin  Kessler  111,  Geophysics  Research 
Directorate,  October  1957. 

No.  57.  Mean  Monthly  :\0(>  mb  and  200  mb  Contours  and  500  mb,  300  mb,  200  mb  Temperatures  (or  the 
Northern  Hemisphere,  K.  W.  Wuhl,  April  1958. 

No,  58.  Theory  of  Large-Scale  Atmospheric  Diffusion  and  its  Application  to  Air  Trajectories,  S.  li. 
Solot  and  K.  M.  Darling,  Jr.,  June  1958. 
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